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ABSTRACT 


The newsboy problem solutions under the conditions of 
risk and of uncertainty about demand are well known. The 
former is the case where the distribution of demand is known 
or estimated, and the latter is the case where it is not 
known but the range of demand is given. 

Minimizing the expected cost under risk and the minimax 
approach under uncertainty are well-known methods to solve 
the problem. The situation is considered where demand fre- 
quency data 1S acquired, one observation per decision period. 
For the expected value solution, this study presents five 
Candidate estimators based on the order statistics and 
evaluates them by comparing costs with those achieved using 
the minimax rule under uncertainty. This is done by simu- 
lation; the results also provide information about the switch- 
ing period from the minimax rule to the proposed quantile 


estimators. 
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I. INTRODUCTION 


The newsboy problem was described by Morse and Kimball 
early in 1950 [1]. This problem has been presented in the 
literature under a variety of names, including newsvendor 
problem, christmas tree problem, and high-fashion apparel 
problem. This kind of problem may be described in the 
following way: The newsboy has to decide how many papers to 
purchase for resale. If he buys too many papers he will 
incur loss due to left-over papers (little salvage cost). 
If, on the other hand, he buys too few, he should also incur 
an opportunity cost due to lost sales. To determine the 
Optimal order quantity, it is necessary to balance two 
Opposing costs. Solutions to the newsboy problem under both 
conditions of risk and uncertainty are well known. The 
determination or estimation of the demand distribution is 
Critical in dealing with this problem. Optimal solutions to 
the newsboy problem which minimize expected cost require 
information about a particular quantile of the demand dis- 
tributions. If the demand follows some known distribution 
the problem can be solved by decision making under risk, as 
is also the case when the key quantile is known. With only 
information concerning the range of demand, optimal solutions 
can be found by some procedure such as the minimax rule under 
uncertainty. 

The purpose of this study may be summarized as follows. 


We are concerned with decision policies during a period when 





demand data is being acquired at the rate of one observation 


a eet ce rT 


per decision period, starting with no data except an estimate 
of the maximum demand. To find the optimal expected value 
decision rule, five candidate quantile estimators will be 
proposed on the basis of order statistics. They will be 
evaluated by comparing their performance with each other and 
with the well-known minimax rule under uncertainty. An esti- 
mator which has the smallest relative cost deviation to the 
risk cost solution (ideal case) will be proposed as optimal. 
We shall also examine the switching time to change from the 
uncertainty rule to applying the proposed estimator. 

In Chapter II, we shall review the newsboy problem solu- 
tions under conditions of both risk and uncertainty, and then 
Suggest alternative approaches for demand estimation, and 
for the switching period from uncertainty to risk. For esti- 
mation, five quantile estimators will be proposed in Chapter 
III. In Chapter IV, candidate quantile estimators will be 
evaluated by simulation. The performance of the estimators, 
and candidate switching periods will be discussed in Chapter 
V. Conclusions and recommendations, for further work, will be 


given in the final chapter. 





Pit oreo OF THE PROBLEM 


In this section we shall briefly describe the well-known 
newsboy problem, illustrating its optimal solutions under 
both risk and uncertainty with respect to demand. Then we 
shall describe the quantile estimation problem we wish to 


solve. 


A. THE NEWSBOY PROBLEM 

The problem may be generalized in the following way. 
The decision maker selects a number (the newsboy's order 
quantity) but the actual demand (the number of papers his 
customers will want to purchase) may have a randomly deter- 
mined value. If the decision maker's number is greater than 
this he pays a cost proportional to the difference, or he 
pays a cost proportional to the amount by which his number is 
smaller than this. His cost is zero only if his number is 
the same as the actual demand. 

Let demand D be a continuous random variable and let 
moantity S be the number the decision maker selects to have 
on hand. The newsboy cost equation may be formulated as a 


two-piece continuous linear function in the following manner. 


C(S - D) one 0 


[A 
O 
[A 
— 


C(S) = CE) 


C,(D- S$) set Dass >, 





where C. is unit cost of surplus when S exceeds D, and Co 
is unit cost of outage or shortage when D exceeds S. The 


cost function is shown in Figure l. 


€ost €(S)} 





0 S Demand, D 


Figure 1. The Newsboy Cost Function 


feet, NEWSBOY PROBLEM SOLUTIONS UNDER RISK 

For our case where demand is a continuous random varia- 
ple D with probability density function f(D), and § is a 
decision variable whose value the decision maker is trying 


to find over D, the expected cost may be represented as: 


S co 
meteiee= C,) (S-b)t(p)dp +c }{ (D-s)f(D)dD- 
0 S 


(2) 
We wish to find quantity S, which minimizes the expected 
cost. The optimum (minimum in this case) value of S may 
be found by differentiating the expected cost function (2) 
With respect to S, setting the derivative to zero, and solving 


for the optimum value S?*: 





00 set 


dE(C(s)] _ : : = 
| ne ) 2a ee i £(D) dD = 0. (3) 


From this we obtain the optimal value S* using the cumula- 


meevye distribution function 


x 
Reem = i 2(D) dD, 
0 


and the result is the well-known result, 
rr (4) 


Differentiating Equation (3) with respect to S to get the 


second derivative at S*, 


aE [C(S) ] 


eeG iC ef (5%) 
ds? | sass email 


Hence, S* is the minimum expected cost solution since 
(C_+C.)£(S*) 1s positive. 

To summarize, the optimal solution S* to minimize the 
expected cost is the [Cf (Cl, +C,) Ith quantile of the demand 
distribution, as shown in Figure 2, and the crucial infor- 


mation we need about the distribution is that quantile. 


C. THE NEWSBOY PROBLEM SOLUTIONS UNDER UNCERTAINTY 
The decision maker may be confronted with a newsboy 


problem without knowing or being able to estimate the 
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0 or Demand, D 


Figure 2. Minimum Expected Cost Solution 


probability distribution of demand. The problem of deciding 
on S in this situation is called a decision under uncertainty. 
The decision maker may be able to estimate an upper bound 
for the demand value. In this case we say that all we know 
about demand is that demand D is such that 0 <« D< Dy. 
“. Three well known approaches to decision making under 
uncertainty are Laplace solution under assumption of a uni- 
form distribution of demand, minimax cost solutions and 
“minimax regret solutions [2]. For the newsboy problem, 
these all lead to the same rule to find the optimal solution. 
We shall discuss only the minimax cost approach. This 


approach is to choose S so that the worst possible cost will 


be minimized, that is, to: minimize maximum cost. The maxi- 


Mum cost occurs at D= 0 or at D= Dae depending on the 
cost coefficients. This can be seen from Figure l. If 
maximum cost occurs at D = 0, we would want to reduce § 
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wae, 1fL Maximum cost occurs at D= Dax! we would want 


to increase S to reduce cost. We can minimize the maximum 


cost when S is chosen so that the cost at D = 0 is equal to 
the cost at D = Dax’ Equating these, we obtain: 
CS 4a Ca D oS 
Ss O° max 
Gis 
Co 
ke Ree eh eS ; 
S (a +e ) Dae 3 a (5) 
Ss oO 
So far, we have introduced the minimum expected cost 
solution under both risk and uncertainty conditions. The 


minimax approach is tied to our initial estimate of D ; 


ax 


and thus the optimal value S* from the minimax rule will 
be in error to the extent that our estimation of D ax is 
in error. This error will continue over many decision 
periods unless our estimate of D aay 1s revised. 
D. QUANTILE ESTIMATION AND A DECISION RULE SWITCHING 

POINT FOR THE NEWSBOY PROBLEM 

We are interested in finding a quantile estimator which 
gives better results than the minimax approach under condi- 
tions of limited demand data, and also in finding the 
Switching period which is the time at which we switch from 


the minimax rule to the expected cost rule using the quantile 


estimator. The purposes of this study are to find the “best” 


ea 





quantile estimator among five candidate estimators, and to 
obtain information about the switching period. 

We recall that we wish to find the quantile of the 
@umutative distribution function, F(D). The optimal S* is 
found by an inverse mapping of the [C./(C,+C.) I] th quantile 
of F(D) into demand D under risk condition. Let q = Co/ 


(C +C_), then we have: 
S Oo 


Although we don't know the demand distribution, each suc- 
cessive decision period provides another observation from the 
demand distribution, and we will assume that an observation 
of demand is indeed obtained each period. As we acquire more 
data we might be able to estimate S* = Eten as a function 
of the given data. Two kinds of approachs are: 

1. a parameteric approach to hypothesizing the form of 


Wn 


demand distribution, 


feed Non-parametric approach to find the quantile. 

For the first case, the normal, and exponential distribu- 
tions have been found to be of considerable value in des- 
cribing demand distributions. The normal distribution has 
been found to describe many demand functions at the factory 
level; the exponential, at the wholesale and retail levels 
[3]. Of course, these distributions should not be auto- 


matically applied to any demand distribution. Statistical 


LS 





tests should establish the basis for any standard assumption 
concerning a demand function. 

feeere directing our attention to the case where the 
expected value solution is not applicable because of lack of 
available information. Further, if we had enough information 
for a fit of a distribution, we should have a good estimate 
of the quantile we seek. Hence, the second case, the non- 
parametric approach, will be considered to find the optimal 


pemwana switching period. 
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III. DEVELOPMENT OF THE ESTIMATORS ON 
THE BASIS OF ORDER STATISTICS 

As mentioned in the previous chapter, we are subject 
to find the optimum (minimum cost solution) of a newsboy 
problem without knowing the q-th quantile of the demand dis- 
tribution. Consequently a non-parametric or distribution-— 
free approach is suggested to find the solution. In this 
chapter, we will introduce the notion of order statistics as 
a basis for estimating this quantile, and suggest five can- 


didate quantile estimators which may be used in newsboy 


Solutions. 


Peto LlMATORS USING ORDER STATISTICS 

The non-parametric approach we shall follow is on the 
basis of order statistics and although it is applicalbe to 
both continuous and discrete random variables, we shall direct 
our attention to the continuous case. We are interested in 
the case where the population is not known, but we have a 
random sample of size n from its unknown cumulative distribu- 
meen function F(X). Let xX,,X 


ee Die  ! 
Meteee 4) < X(2) < ... < X{n), where the X(1) are arranged 


xo be a random sample. 


in order of increasing magnitudes, and are defined to be the 
Order statistics. Let Be be the G=th quantile” of the popu- 
Jation. The estimate of a can be expressed as a function of 


Ky rXoreee Xe ie. ; 
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Estimate of a = £(X),X x) 


nrees7 


We shall consider now a function F(x) which is construc- 


prxgre ce kX: The function 


Fx) 1s called empirical cumulative distribution function. 


ted from the sample values of X 


Let Kp rXore eer KX, denote the observed values of Xp rXoree eX 
Then for each number x(-~ < x < ~), the value of F(x) is 
defined to be the proportion of observed values in the sample 
Bitch are less than or equal to x. It follows from the 
Glivenko-Cantelli Theorem that, as n approaches infinity the 
empirical cumulative distribution function Fs) converges 
uniformly to F(x). Since we are concerned with an unknown 
cumulative distribution function F(x) which is continuous, 
some type of smoothed curve of Ff) might yield a reasonable 
estimator of F(x). Let q be the [Cc /(C, +O.) Ith quantile, 
where we wish to find the estimate of _ so that 

F(estimate of SE = q. This is represented in Figure 3. 

Now, we are interested in the estimator which gives us mini- 
mum cost, but we have no infcrmation about it yet. Let us 
turn our attention to the expected cost as the function of the 
estimate of ae Let perturbation be the difference between 
the estimate of Se and true value of x Then the expected 
Peet Of perturbation is proportional to the variance of the 


given unbiased estimator [4]. The properties of unbiasedness 
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X-1) Xi) X tt) Demand, X 


Figure 3. Estimate of X_ on the Smoothed 
Curve of Empi¥ical Cry, FL tx) 


and minimum variance should be considered for selecting 


estimators. 


B. RATIONALE FOR THE ESTIMATORS 

In this section, we shall suggest five different quantile 
estimators on the basis of order statistics. Ina statisti- 
cal evaluation of these estimators, it is desirable to use 
the uniform (0,1) distribution, since the expectation and 
variance of its order statistics are simple and exact. 
Investigating these properties of order statistics using other 
distributions is not easy. If we have a distribution similar 
to the uniform, or if we have a large number of observations, 
some approximations may be possible. Three points which are 
near each other can be assumed to be linear and some estimator 


using order statistics would be approximately unbiased. 


ley 





We shall examine the properties of the order statistics 
from the uniform (0,1) distribution, and we shall consider 


three neighboring order statistics. Let X Xerys and 


a=.) 


Xr t1) be the (r-l)th, r-th, and (r+l)th order statistics 


respectively. Define the value r as: 
r = [ng + 0.5] , (6) 


where [X] denotes the largestinteger of X. Let us consider 


an estimate of a as the linear combination of those; 


estimate of X = aX + bX + cX 

( q? (ee) (r) (r+1)’ 
where a, b, and c are nonnegative and sum to unity. There 
can be many possible ways to take estimators, but we shall 
examine typical ones considering the unbiasedness and the 
magnitude of variance. 


1. The First Candidate Estimator: a 


First of all, the estimator should be unbiased, 


ree. , Elles = a If we take the r-th order statistic, 
ee) as the first candidate estimator, 1.e., 
Mee 7 


this is an unbiased estimator since E[X iv—=er/ (ntl) for 


iat) 


the r-th order statistic. The coefficients are a = c = O and 


io 





b= 1in this case. The first candidate estimator is repre- 


sented in Figure 4. 


F (x) 





‘ 
~~ ween ee &e Se S&S | ee = 
‘ 


Demand, X 


*(r-1) Kr) r+) 


~ 


Figure 4. The First Candidate Estimator, - 


~ 
~ 


g 
We shall propose (r-1)th and (r+l)th order statistics 


2, The Second and Third Candidate Estimators: Xo 
as the second and third candidate estimators: 
x ra Ririt) ’ (8) 


and 


bd te 


g = KR r41)- (9) 


Even though these are not unbiased estimators, we are not 
certain which one gives us better performance since we are 
considering the properties on the uniform distribution and 


there is some round up error for taking the value r. 
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3. The Fourth Candidate Estimator: Xo 
We shall propose an equally weighted linear combina- 
Eronmot the (r-1)th and (r+l)th order statistics as the 


fourth candidate estimator: 


1 
x lx 


| 


g ae + Xray) I> C0) 


This estimator is an unbiased estimator, since 


and has the minimum variance using 


r(n-s-l) 


COV wie |= 5 
(n+1)° (n+2) 


(r) (s) 
of the uniform distribution [5]. The variance of this esti- 


mator is given as follows: 


+ *(x41), Zz (n-r) (2r-1)+(r-1L) 


areas 
Var[x ] = var (122 Le 5 
q 2(n+1)* (n+2) 


This estimator is represented in Figure 5. 
4. The Fifth Candidate Estimator: Ke 
We shall propose an equally weighted linear combina- 
tion of the (r-1l)th, r-th, and (rt+l)th order statistics as 


the fifth candidate estimator, or 
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Demand, X 


(cr) *(r41) 


Figure 5. The Fourth Candidate Estimator, Xy 


ae: 
Clee 3 | 


>< | 


+ X aa (11) 


esi Pe Get) 


This is also an unbiased estimator and the variance of this 


estimator is given by: 


var [X_] = Anor) (rns) terns) 
i 9 (ntl) * (nt2) 


The difference in variance between the fourth estimator and 


in = 
T8(n+1l) (nt2)° When n = 1 and 
4 


Peewl0) this difference is approximately 0.09 and 4x10 ", 


the fifth estimator is equal to 


respectively. The difference converges rapidly to zero as n 
increases. Since the quantity difference may be assumed to 
be negligible and the variance is only concerned with uni- 
form distribution, it is hard to say which estimator gives 


us better performance. 


ZL 





a 


In summary, although there may be many possible 
candidate estimators depending upon different forms of 
linear combination, it may be reasonable to examine five 
previously proposed estimators. They can be shown with 


properties as: 


an 


ieethe First Estimator: a = X ery (unbiased) 
2) The Second Estimator: xX = xX, (biased) 
q (a= 1) 

3) The Third Estimator: a = Set) (biased) 

he —_ 1 | ie: 
4) The Fourth Estimator: a = RIK gy * Xray} 

(unbiased min. variance) 
: ; = 1 
: = = + + 

meee Fitth Estimator Ae z1X (p17) Rr) Siaye 


(unbiased). 


Having proposed these quantile estimators, it is of 
interest to see how well they perform in the context of 
newsboy decision making, under cost criteria. This is the 


subject of the next chapter. 


22 





IV. EVALUATION OF THE CANDIDATE ESTIMATORS 


In this chapter, the five candidate estimators will be 
tested and investigated for their performance in the newsboy 
problem for several demand distributions. As performance 
measures, two kinds of costs, namely an average cumulative 
cost and a mixed cost which we shall introduce later on, will 
be obtained during the simulation and then the performance 
of the candidate estimator will be evaluated using the rela- 


tive cost deviation to the ideal cost. 


A. DESIGN FOR SIMULATION 

A newsboy simulation will consist of 50 replications, 
each of which simulates 50 decision periods for a given 
demand distribution. The quantile value (q = Co/ (C5 +C0)) 
for each distribution is varied from 0.1 to 0.9 in increments 
CreanO. 2. 

To find the value r and consequently to find the estimate 
Ge Kae some modification of the method of determining r is 
necessary during the first and the last few periods, since 
there is an insufficient amount of data to apply the candi- 
date ec imators. For example, if the value r for finding 
the estimate is equal to zero, then Xo1) will be chosen as 
the estimate of Xa and if r =n, then the estimate of xX, 
will be Xen) and so forth. 

For each period, the average cumulative cost and the 


mixed cost are computed for seven different cases which are 
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represented as follows: 
1) Risk case (Ideal solution), 
2) Uncertainty case using the minimax rule, 
3) <A case uSing the estimator XQ, 
4) A case using the estimator Kae 
5) A case using the estimator ko, 


6) A case using the estimator Xo and 


7) A case uSing the estimator Xo 


The cost for each period is averaged over 50 replications, and 
then the average cumulative cost at period k is computed by 
cumulating the average cost from the first period to k-th 
period. The mixed cost at k-th period is given by summing 

the average cost using the minimax rule up to (k-1)th period 
and the average cost using an estimator from the k-th period 


to the last period. 


B. DISTRIBUTION AND PARAMETERS 

As mentioned previously, one source has suggested that 
the normal distribution may be considered as the demand dis- 
tribution at the factory level and the exponential distribu- 
tion may be considered as the demand distribution at the 
wholesale and retail levels, and so on. Five characterized 


forms of demand distributions including symmetric, right-skewed, 
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and left-skewed will be selected to test the performance of 
the candidate estimators for generality, and are shown in 
Figure 6. 


£7Cs¢) 


eyo NO 
LP ae 


0 Demand, X 


Figure 6. Characterized Forms of Demand Distributions 


Demand distributions used in the simulation and their 
parameters are shown in Table I. The parameters for the 
normal distribution were selected so that negative demand 


Should not be generated. 


Table I. Distributions and Their Parameters 









Zo) 





For each demand distribution, the 99.9th percentile of 
demand except uniform distribution will be adopted as Da 
to apply the minimax rule under uncertainty. Quantile values 


g = Co/ (CS + CS) wise 050.3, 0.5, 0.7, and 0.9 for each 


@estribution. 


C. ALGORITHM FOR SIMULATION 
The program steps are ordered in the following way: 
1. Define a new demand distribution. 
2. Assign the quantile value and maximum demand. 
3. Compute the optimum in cases of uncertainty and risk. 
4. Generate three random deviates from the specified 
demand distribution, representing demand for the first 
three periods. 
5. Obtain the order statistics. 
6. Compute the estimates using the candidate estimators. 
7. Generate a random deviate for the specified demand 
distribution for the next period, compute the cost for 
the period and save it. 
8. If the period is beyond the final period, go to Step 9; 
otherwise go to Step 5. 
9. If the number of replications is greater than the maxi- 
mum, go to Step 10; Otherwise go to Step 4. 
10. Find the average cost for each period over total 
replications. 
11. Compute the average cumulative cost, the mixed cost, 
and their relative cost deviation to the risk cost. 


Then tabulate the results. 
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ieeebt quantile value is greater than 0.9, go to Step 13; 
otherwise go to Step 2. 
ieeeLe we want to try another demand distribution go to 


Step 1; otherwise stop. 


D. DISCUSSION OF THE SIMULATION RESULTS 

In this section, the simulation results will be displayed 
and evaluated. As mentioned previously, the simulation run 
was performed over 50 replications, each of which consisted 
of 50 consecutive demand periods. From these trials, the 
average cumulative costs and the mixed costs were computed 
for various estimators. 

Table II shows the average cumulative cost for 50 periods 
for the given quantile values and specified demand distribu- 
tion. Errors in estimating Dae were not considered. This 
performance of quantile estimation may be compared directly 
with "best case" minimax results. Within a distribution we 
IMay conclude one estimator is better than the other but it 
is not easy to conclude for all distributions using Table II. 

Now, a measure of effectiveness (MOE) to evaluate the 
performance of the estimators for all distributions is needed. 


Let us define the relative cost deviation of an estimator as 


follows: 
Belative cost deviation _ (cost using estimator) - (risk cost) 
(to the risk cost) (risk cost) 


Ze 





tab lew f 


Average Cumulative Costs for 50 
Periods and 50 Replications 
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aes c * ni Estimator 










Quantile (q) 
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20) X 
~q 
xX 
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4Q 





eae GHZ Zleleal Zee cn 12.496 






Risk 6 
Gamma Minimax 7.024 23.094 36.970 37.240 18.098 
; - Bit. | 7.253 16.963 22.275 22.141 13.649 
| x is Cones OG 2442) de iise 
i 7B U7 jEI~ DAG Bea acy We cey 
x" 7.305 17.067 22.100 21.992 13.863 
Lz ao OA Oc ene G 4 s.299 
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Table II (Continued) 


Desc Quantile (q) 
bist i Ow 0.5 


Sil ><] Poot 
Oo QM AQ 


1Q 
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This relative cost deviation may be used as an MOE to evalu- 
ate the performance of an estimator. The best performance 
was given by the risk solution during the simulation, an 
expected result with regard to theoretical considerations. 
Table III exhibits the relative cost deviation for 
quantile values from 0.1 to 0.9 for each demand distribution 
investigated. As a special case, if we pursue the uniform 
demand distribution in Tables II and III, the risk and 
Minimax approach would match each other because the risk 
cost solution Ko Such that F(X) = Co/ (C5 + S) is the same 
as the minimax cost solution given by [Cc /(C,+C,) 1D, 


Gan 


whenever the maximum demand dD. 1s determined exactly. If 


ax 
we read through the third column of overall distribution in 
Table III, an estimator, Xa shows the best performance among 
candidate estimators at quantile 0.1. Also, Xa did better 
at quantile 0.9. Excluding these two points, the estimator 
Xa gave the best performance at all other quantiles. Within 
the uniform distribution the minimax approach shows the best 
performance (since it is identical to the risk solution), 
but it shows the worst performance in the other distribution. 
Similarly, we can say which estimator represents the best 
performance within a distribution by reading through row 
by row, and which estimator shows the best performance at 
a quantile by reading through column by column. 

We are interested in the estimator which gives us the 


best performance for all quantiles over all distributions. 
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Table. © bal 


Relative Cost Deviation From The Average Cumulative Cost 


i 
Dist'n 
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Risk 0.0000 0.0000 





0.0000 

































































Uniform |Minimax | 0.0000 0.0000 0.0000 
a) x, 0.0509 0.0696 0.5104 
X 0.0634 0.1066 0.9909 
x 0.0517 0.0496 0.5844 
% 0.0318 0.0493 0.5655 
x 0.0345 0.0502 0.4960 
Risk 0.0000 0.0000 0.0000 0.0000 0.0000 
Normal 0.6149 0.4229 0.0405 0.1154 1.5526 
(35,100) 
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020877 
Q20558 
0.0459 
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Oe S23 
NeSoo2 
OEZOL 
0.1307 


0.0000 
Peso 
0.6388 
1.2543 
one 
uz 
Gato 


0.0000 
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2.1542 
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S24738 
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One way to compare is to note that the performance of an 
estimator overall quantiles over all distributions can be 
represented as the total sum of the relative cost deviations. 
Looking at the last column for overall quantiles over all 


distributions it may be concluded that the best performance 


SE 
Been) © (x) 


+ yee) | for all quantiles over all distributions as shown 


(minimum sum) is given by the estimator Xa = 


in Table III, because the minimum relative cost deviation is 
the closest to the risk cost (ideal cost solution). For 
specific quantiles, of course, occasionally other estimators 
performed better. In an application of the newsboy problem 
under initial uncertainty, one has information about the 
quantile q but often not about the form of the distribution. 
meearaingly, the information at the foot of Table III may 


provide information in selecting an estimator. 


E. DISCUSSION OF SWITCHING RULES 

When solving a newsboy problem repeatedly under uncer- 
tainty conditions we should probably use the minimax rule for 
the first period since there is no applicable amount of data 
to use in quantile estimation. Each successive decision 
period provides an additional piece of demand data which 
Will contribute to quantile estimation. Therefore we are 
interested in the number of decision periods that should 
elapse before we apply the proposed estimator Xo in place of 
the minimax rule. From the simulation results, Table IV 


shows the mixed cost with switching at period 25, the middle 


oS 





Table IV 


The Mixed Cost for the 25-th Period 
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Table IV (Continued) 


Desc. 
Sa Estimator 


x 
4 
x 
=7 
X 
ao 
Xx 
<A 
X 
q 
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period in our fifty-period simulation. The mixed costs 
here are the sum of costs for the first 24 periods given by 
the minimax rule and the costs given by an estimator from 
the 25-th to the last period. When we examine various dis- 
tributions in Table IV the minimum mixed cost is obtained 
by the risk approach, which is to be expected, and the 
proposed estimator, Ko appears to provide the minimum 
mixed cost among candidate estimators, although by a narrow 
margin. 

Now, we wish to find the number of periods to retain 
the minimax rule before applying the proposed estimator Ko 
Mixed costs were computed for various switching periods to 
find the period which gives the minimum mixed cost for each 
quantile of every distribution except uniform distribution, 
Since the proposed estimator Xe doesn't give better results 
than the minimax rule for the uniform distribution. Table 
V shows the best switching periods (minimax to estimator X,,) 


for the quantiles and distributions used in the simulation 


during 40 periods. 


Table V 


Best Switching Periods Using the Proposed Estimator Xa 


Quantile (q) 
O23 Oy 


Normal 


Exponential 


Gamma 


Beta 
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Looking at Table V, we see that the best switching period 
momez7eat Guantile 0.5 of the normal distribution. This 
result is caused by the fact that the minimax rule behaves 
well around the median when demand has a symmetric distri- 
bution. The occurrence of switching period 48 at quantile 
0.9 of the beta distribution can be explained by the fact 
that with the beta distribution used in the simulation, 
the optimal S* obtained by the minimax rule is close to 
the ideal S* at this quantile, i.e., the minimax S* at 0.9 
imoemecual to 0.8767 and the ideal S* is 0.8660. 

We shall suggest 6 periods, 4, 5, 6, 7, 1l, and 13 as 
the candidates for the switching period for all quantiles 
over all distributions. Later periods are not included in 
the candidates since the estimator Ko shows much better 
performance (with reference to the foot of Table III for 
nearly all quantiles over all distributions). Table IV 
shows the mixed costs obtained by the estimator Ka for the 
25th period. Similar tables (not shown) were computed for 
various other periods, and showed that the mixed costs 
increased as the period increased. We propose the relative 
mixed cost deviation for a candidate period using the esti- 
mator x as an MOE. It may be reasonable to conclude that 
the period giving the minimum is the most preferable. Table 
VI shows the relative mixed cost deviation at each candidate 
period for all quantiles over all distributions using the 


estimator Xo 
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Table VI 


Relative Mixed Cost Deviation Using the Estimator Xo 






Switching 
Period Distribution] Estimator 






Quantile (q) 
O25 On 0.9 Total 
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5 Exponential .0548 .0851 .3337 
0422 .0708 .2971 
.1003 .0868 .4578 
fc 5606 1.9793 
.0457  .0626 .0693 
.0863 .0663 .0644 
6 .0697 .0742 .0765 
0427 0675 .0814 
.1591 .1400 .0790 
.4035 .4106 .3706 
Uniform .0300 .0420 .0547 .0601 .1020 .2888 | 
Normal | 1311 .0846 .0533 .0572 .1084 .4346 | 
7 | Exponential | x, .0738 .0924 .0980 .0944 .0987 .4573 


Gamma | .0423 .0487 .0754 .0924 .0847 .3435 
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Table VI (Continued) 


Switching 
Period Distribution | Estimator Quantile (q) 
Gas G5 O27 
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The relative mixed cost deviation at each candidate 
period is shown in Table VII. The minimum deviation is 
given at period five and the deviation is increasing as the 
period increases from the fifth period. Using this result 
we May propose that the switching period from the minimax 


rule to applying the estimator Ka Should be the fifth period. 


Table VII 


The Relative Mixed Cost Deviation 
At Each Candidate Period 


2.0346 





1.9793 (minimum) 





Zao ge2 







2.3174 


3527354 






3.9654 
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V. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 


The minimax approach is a well known approach to decision 
making under uncertainty. For a newsboy-type problem, we 
developed five candidate estimators which could be used as 
data aggregates, and evaluated them by comparing with the 
Minimax rule. These candidate estimators were tested for 
aiemerles, 0.1, 0.3, 0.5, 0.7, and 0.9 of several distribu- 
tions: uniform, normal (symmetric), exponential, gamma 
(right-skewed) and beta (left-skewed). Since all possible 
quantiles and other kinds of distributions were not covered, 
there may be some amount of loss of generality. 

With these limitations, as discussed in Chapter IV, the 
mit) 8 ie) FSi) developed on the basis 


of order statistics showed the best performance among those 


estimator Xo = 


considered. The simulation results suggest that if a single 
decision procedure over all quantiles is needed, the esti- 
mator Ka is recommended and the time for switching to this 
estimator should be at the fifth period. The proposed 
meeeeaure may be summarized in the following way: 

(1) Estimate Dax’ 

(2) Compute the guantile gq = Co (C5 + C,)- 

(3) Apply the minimax rule S* = Cae POE Che. fairs t 

four periods. 


(4) Compute the value r such that r = [ng+0.5] from the 


fifth period where n is the period at which we wish 


41 





to estimate the demand and [X] denotes the largest 
integer of X. 

ory) 98S +1) 
and apply the proposed estimator XQ 


fy Compute the order statistics X 


II 
i 


a 
q Bee eee 


In this study, it was implicitly assumed that the value 
of demand for each period was obtainable even if D> S. For 


further study, the case when demand data greater than the 


quantity S is not available is worthy of investigation. 


pennant 


—— 





Also, a non-parametric approach which can handle extreme 


quantiles (e.g., 0.01 or 0.99) and extremely skewed distri- 
Mons could be my eet The comparison between a 
parametric approach when the demand data are assumed to fit 
a probability distribution and non-parametric approach is 
also recommended. Finally, future simulation work on this 


tical analysis could be made of the simulation results. 
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